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The mass-metallicity relation of interacting galaxies 
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ABSTRACT 

We study the mass-metallicity relation of galaxies in pairs and in isolation taken from 
the SDSS-DR4 using the stellar masses and oxygen abundances derived by Tremonti et al. 
(2004). Close galaxy pairs, defined by projected separation r p < 25 kpc h^ 1 and radial ve- 
locity AV < 350 kms -1 , are morphologically classified according to the strength of the 
interaction signs. We find that only for pairs showing signs of strong interactions, the mass- 
metallicity relation differs significantly from that of galaxies in isolation. In such pairs, the 
mean gas-phase oxygen abundances of galaxies with low stellar masses (M* < 10 9 M.Qh~ 1 ) 

exhibit an excess of 0.2 dex. Conversely, at larger masses (M* > 10 10 M Q /i _1 ) galaxies have 
a systematically lower metallicity, although with a smaller difference (—0.05 dex). Similar 
trends are obtained if g-band magnitudes are used instead of stellar masses. In minor inter- 
actions, we find that the less massive member is systematically enriched, while a galaxy in 
interaction with a comparable stellar mass companion shows a metallicity decrement with 
respect to galaxies in isolation. 

We argue that metal-rich starbursts triggered by a more massive component, and inflows 
of low metallicity gas induced by comparable or less massive companion galaxies, provide a 
natural scenario to explain our findings. 

Key words: galaxies: formation - galaxies: evolution - galaxies: abundances - galaxies: in- 
teractions. 



1 INTRODUCTION 

Chemical features observed in galaxies can store fossil records 
of their history of formation since they are the result of differ- 
ent physical mechanisms acting a t different stages of evolution 
dFreeman & Blan d-Hawthorn 2002). Among them, galaxy interac- 
tions and mergers are considered important processes which can af- 
fect the star formation activity significantly as has been reported by 
numerous works (e.g. Sersic & Pastoriza 1967; Tinslev & Larsonl 
19781: b arton et aj] l200Cl : lLambas et al"]|2003l ; iNikolic et alj|2004 



Alonso et alj 200t ). The link of interactions and mergers with 



chemical properties have been recently started t o be addressed 
(iDonzelli & Pastorizal |2000| : iMarquez et alj |2002| ; I Fabbian o et al.1 
l2004lKewlev et al.l2006h . Recentlv. lKewlev et alJd2006t) found the 
metallicity in the nuclear regions of galaxies in very close pairs to 
be displaced to lower levels than those measured in more distant 
pairs and isolated galaxies. The authors claimed this effect to be 
caused by gas inflows triggering by the interactions. 

In this letter, we study the mass-metallicity relation (MZR, 
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iLequeux et al.1 Il979h of close galaxy pairs as a tool to study 
the effects of i nterac tions on metallicity. In the local Universe, 
iTremonti et alj {2004}) have confirmed the dep endence of meta llic- 
ity on stellar mass with high stastistical signal. lErb et al.1 <2006f) has 
extended the study to high redshift finding a simi lar correlation, 
althou gh displaced to lower metallicity (see also iMaiolino et al.1 
2007). The luminosity-metallicity relation (LMR) has been also 
well estimated in local Universe when stellar masses are not pos- 
sible to measure. However, it has been show n that luminosit y is a 
less fundamental parameter than mass (e.g. lErb et alj |2006). De- 
spite these new results, the origin of the MZR and its evolution 
remain unclear. Mergers and interactions could play an important 
role in shaping the MZR. Hence, a statistical analysis of the effects 
of interactions on chemical properties could help to shed light on 
this issue. 

Numerical simulations showed how tidal interactions can 
trigger star forma tion activity and can modify the morphology 



of galaxies (e.g. | Ke nnicutt et al. 



Miho s & Hernquistl 1 1 9961 : iBarnej |2004|) . Gas inflows gener 



19981 a nd references therein; 



ated during these interact i ons feed star f o rmation activi t y (e.£ 
barnes & Hernquistl 1 19961 : iTisseral [2000). Tissera et al.1 ^2005 
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Figure 1. Mass-metallicity (left) and Luminosity-metallicity (right) for galaxies in M (black solid lines), T (green dashed lines), N (red dotted-dashed lines), 
and D (blue long-dashed lines) type pairs, and in the control sample (magenta dotted lines). The thick solid line corresponds to Tremonti et al mean values. 
Error bars have been estimated by applying the bootstrap resampling technique. 



showed that simulated galaxies in the concordance A-CDM cos- 
mology reproduce the shape of the MZR but not its evolution. 
Part of the failure could be caused by th e lack of Supernov a 
(SN) feedback in these sim ulations. In fact. iBrooks et alJ d2007t) . 
iGovernato et al] d2007l) . and lFinlator & Pavel d2007t) . resorting to 
ad-hoc parametrizations for SN feedback model, claimed to be able 
to reproduce the observed MZR in a cosmological context. 

We performed a statistical analysis of the chemical properties 
of galaxies in pairs selected from the SDSS-DR4, focusing on the 
MZR. Since it is already clear that proximity is the main parame- 
ter associated to important increases in the star formation activity, 
hereafter, we will centre our analysis on close galaxy pairs. These 
pairs have been reclassified according to the level of morpholog- 
ical disturbances detected b y inspection of the ir images in merg- 
ing, tidal and non-perturbed dAlonso et al.l2007h . Merging pairs are 
those exhibiting the higher level of star formation activity. 

This letter is organized as follows. In Section 2 we described 
the galaxy samples. Section 3 discussed the results. The conclu- 
sions summarize the main findings. 



2 DATA 

The main SDSS galaxy sample is essentially a magnitude limited 
spectroscopic sample r« m < 17.77 with most of the galaxies span- 
ning a redshift ra nge < z < 0.25 with a median redshift at 0. 1 
( Strau sTet al.l2002f) . Following lAlonso et al]j200r3 . and references 
therein), we selected galaxy pairs by requiring members to have 
a relative projected separation r p < 100 kpc h^ 1 and a relative 
radial velocity AV < 350 km s _1 . 

In order to properly assess the significance of the results ob- 
tained from the pair catalogue, we defined a control sample by a 
Monte Carlo algorithm that selects for each galaxy in a pair, two 
galaxies with similar redshift and r-band luminosity but without a 



close companion within the adopted relative velocity and separa- 
tion thresholds. This procedure assures that the control catalogue 
has the same selection effects than the pair sample, and conse- 
quently, it can be used to estimate the actual difference between 
galax ies in pairs and in isolatio n. 

lAlonso etail feood 120071) discussed possible effects of aper- 
ture and incompleteness, which were found to have no significant 
consequences in the analysis. Active galactic nuclei (AGN) could 
contribute to the emission spectral features affecting our interpreta- 
tion. Hence, AGNs had not been included in neither the control nor 
the pair catalogues. 

As mentioned above, we focus our analysis on the close 
galaxy pairs defined as those with a relative projected separa- 
tion r p < 25 kpc h^ 1 . These pairs are subclassified accord- 
ingly to the level of morphological disturbances in non-disturbed 
(N -type), tidal (T-type) and merging (M-type) pairs as discussed 
by lAlonso et al. d2007t) . M-type pairs show clear and strong sig- 
nals of interactions and have active star formation. For complete- 
ness, we have also included the estimations for more distant pairs, 
50 kpc h,- 1 <r p < 100 k pc h' 1 (D-type) 

For the SSDS-DR4. ITremonti etal] d2004l) estimated gas- 
phase oxygen abundance, 12+log O/H, and derived stellar masses 
(M*). In Table Q] we show the numbers of pairs and galaxies in 
the different samples which have oxygen abundances measured for 
both members or for only one member. Pairs with only one member 
with oxygen estimations have been used to improve the statistics, 
when possible. 



3 ANALYSIS AND RESULTS 

In Fig. [T] we show the mean MZR and LZR for galaxies in M, T, 
N and D-type pairs and the control sa mple. We also plot the fit to 
the observed MZR given in Eq. 3 by ITremonti et al] d2004l) . and 
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Table 1. Number of pairs and galaxies in the different samples used in this 
work (N p total number of pairs; Np number of pairs with oxygen abun- 
dances measured for both members; N g number of galaxies with oxygen 
abundances estimated for at least one member) 



Sample 


M 


T 


N 


D 


C 


N p 


383 


686 


536 


4982 




p 


94 


142 


124 


679 




N g 


284 


454 


410 


1358 


6034 



the m ean value of the observed LZR (also given bv lTremonti et al.l 
2004). It can be appreciated that our different samples are in gen- 
eral agreement with these mean relations except for M-type pairs. 
For this sample, at low stellar mass, or low luminosity, there is a 
significant excess in the gas-phase metalli city with respect to th e 
other samples and mean observed values o f lTremonti et all J2004I) . 
Conversely, at high mass (or high luminosity), we find M-type pairs 
galaxies to have slightly smaller met allicity. This behavio ur is in 
agreement with the results reported by Ke wlev et al.l {2006). These 
authors found a decrease in the metallicity in the nuclear region of 
galaxies in close pairs which had members with luminosities con- 
sistent with our high luminosity end. 

We have also explored the dependence of the MZR on galaxy 
concentration for the samples analised finding a similar behaviour 
irrespective of galaxy concentration (here concentration is defined 
as the ratio of Petrosian 90-50 per cent r-band light radii). Since 
the conc entration paramet er is well correlated to morphology (see 
however lEllis et all d2005l) for a detailled discussion of the relation 
of morphology and other galaxy properties), this result indicates 
that our findings are not merely a morphology-induced effect. 

The change in the slope of the MZR of M-type pairs as a func- 
tion of M* shown in Fig. [T] suggests a possible dependence of the 
effects of interactions on the stellar mass of the members. Hence, 
we splitted the pair samples into minor and major interaction sub- 
samples by adopting a ratio of 0.20 between the corresponding stel- 
lar masses. We studied the MZR and star formation activity for each 
component in major and minor interactions in the four defined pair 
samples. 

As shown in Fig. [2] two cases depart co nsiderably from the 
global relation given by Tremonti et al.l d2004l) . Galaxies in major 
interactions (Fig.[2^), which show a lower metallicity, and the small 
component of minor interactions (Fig. |2p) with strong disturbances 
(M-type) which exhibit a significantly larger mean metallicity at a 
given stellar mass. 

These results suggest that in major interactions, there is an 
efficient inflow of pristine gas from the external regions of the halo. 
We argue that the lack of an effect for the large component (M| ) in 
minor interactions may reside on the relatively weak perturbation 
these galaxies are subject to. On the contrary the small component 
(M*) of minor interactions are the ones mostly disturbed by the 
interactions and therefore those subject to major episodes of star 
formation and metallicity enhancement. 

We argue that the ratio between the stellar masses of the galax- 
ies involved in an interaction is an important parameter in setting 
the metallicity content of galaxies. To explore this, we show in 
Fig. [3] (upper panel) that the metallicity difference (A) between 
the less massive and the most massive components in M-type mi- 
nor interactions is significantly reduced in comparison to that of 
galaxies in N-type pairs (similar to D-type pairs). Members of ma- 
jor interactions in M-type pairs have similar metallicity differences 
than those of N-type pairs (lower panel in Fig. [3}- 
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Figure 2. Mass - metallicity relation for the different type pairs: M (black), 
T (green), N (red), and D (blue); in major interactions (M| /M* > 5) 
(a) and in minor interactions (M£/M* < 5) (b: massive components; c: 
less massive components). The thick solid line corresponds to Tremonti et 
al fit. Error bars have been estimated by applying the bootstrap resampling 
technique. 



An analysis of the star formation activity indicates that metal- 
rich less massive members in M-type pairs are in general also those 
exhibiting a high star formation activity, at a given stellar mass. 
This somewhat expected result reinforces the idea that minor com- 
ponents in strong interactions are likely to have experienced strong 
gas shocks leading to star formation and enrichment. 

In order to further assess the role of the relative mass ratio of 
the galaxy members on the mean gas-phase metallicity, we calcu- 
lated the metallicity excess (e) as the difference between the oxygen 
abundance of a given galaxy in a pair and the mean abundance of 
galaxies in the control sample with similar stellar mass as a func- 
tion of the ratio between the stellar mass of the companion (M c ) 
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Figure 3. Distribution of oxygen abundance differences between the less 
massive and the most massive components in M-type (black solid lines) 
and N-type (red dashed lines) pairs in minor (upper panel) and major (lower 
panel) interactions. 

galaxy and the stellar mass of the given galaxy (M g ). We estimated 
e for all subsamples. In Fig. [4] we display the results corresponding 
to pairs in the M-type and D-type subsamples since these galaxies 
set the range of variation. As it can be appreciated, in M-type pairs, 
galaxies with massive companion (log M c /M g > 0) have signif- 
icantly enhanced their metallicity with respect to galaxies in the 
control sample with the same stellar mass. The opposite occurs for 
galaxies with a less massive companion (log M c /M g < 0) where 
the metallicity is lower. Conversely, D-type pairs exhibit a remark- 
able consistency with the control sample. The behaviour of M-type 
pairs stresses that the relative mass ratio plays a major role in deter- 
mining the metallicity properties in galaxies in strong interactions. 



4 DISCUSSION AND CONCLUSIONS 

We analyse the effects of interaction on the metallicity proper- 
ties of galaxies. For this aim we have studied pairs of galaxies 
with projected separation r p < 25 kpc h^ 1 and radial veloc- 
ity AV < 350 km~ a from the SDSS-DR4 and oxygen abun- 
dances from lTremonti et all J20041) . Using the visual classification 
of pairs into M-type, T-type, N-type, and D-type pairs of Alonso 
et al. (2007), we find that pairs undergoing a strong interaction 
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Figure 4. Metallicity excess of a given galaxy with respect to galaxies 
in the control sample with similar stellar mass as a function of mass ratio 
between the companion galaxy and the given one, for M (solid line), and D 
(long-dashed line) type pairs. Error bars correspond to the mean standard 
deviation. 



differ significantly on the mass metallicity relation in compari- 
son to other pair categories and a control sample. In the low mass 
range (M» < lO^o/i" 1 ) the mean metallicity of galaxies in M- 
type pairs is 0.2 dex greater. On the contrary, at large masses 
(M* > lO lo M0/i _1 ), M-type pair galaxies have a systematically 
lower metallicity, although with a smaller difference (—0.05 dex). 

When dividing the pairs into minor (Ml/Mj > 5) and major 
(M2/M* < 5) interactions, we find that the less massive mem- 
ber in strong minor interactions are always enriched irrespective of 
their mass. The interactions of galaxies of comparable stellar mass 
content have a systematically lower metallicity content, suggestive 
of an inflow of unenriched material from the external regions. 

We notice that the finite size of the fiber might affect the re- 
sults in the sense that the smaller the angular size of the objects 
the larger the region where metallicity is derived. Nevertheless, in 
this work we have considered galaxies of similar stellar masses and 
redshifts differing only on the presence of a companion. In order to 
analyse this effect on the MZR into more detail, we have consid- 
ered separately large and small angular size galaxies (Petrosian 90 
per cent r-band light radius ^ 8 arcsec) for the samples in Fig.[T] 
We find that the trends obtained in this work remain unchanged re- 
gardless of the apparent angular size of the galaxies. Therefore, the 
results are robust against a combination of observational aperture 
and metallicity gradient in the galaxies. 

We may interpret our results in terms of a combination of the 
effects of a metal-rich starburst induced by a massive component 
and the pristine gas inflow triggered by comparable or less mas- 
sive companion galaxies. These two effects are expected naturally 
in current galaxy formation and evolution scenarios, so our find- 
ings may help to understand the origin and evolution of the mass- 
metalicity relation. 
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